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Abstract. The signal-to-noisc ratio (SNR) is used in gravitational-wave 
observations as the basic figure of merit for detection confidence and, together with 
Q j the Fisher matrix, for the amount of physical information that can be extracted 

from a detected signal. SNRs are usually computed from a sensitivity curve, which 
describes the gravitational-wave amplitude needed by a monochromatic source of 
given frequency to achieve a threshold SNR. Although the term "sensitivity" is 
s ! used loosely to refer to the detector's noise spectral density, the two quantities 

r- H arc not the same: the sensitivity includes also the frequency- and orientation- 

dependent response of the detector to gravitational waves, and takes into account 
the duration of observation. For interferometric space-based detectors similar 
to LISA, which are sensitive to long-lived signals and have constantly changing 
position and orientation, exact SNRs need to be computed on a source-by- 
> * | source basis. For convenience, most authors prefer to work with sky-averaged 

j^JTJ sensitivities, accepting inaccurate SNRs for individual sources and giving up 

1 1 control over the statistical distribution of SNRs for source populations. In this 

paper, we describe a straightforward end-to-end recipe to compute the non-sky- 
£N| averaged sensitivity of interferometric space-based detectors of any geometry. 

^ This recipe includes the effects of spacecraft motion and of seasonal variations 

■^j- in the partially subtracted confusion foreground from Galactic binaries, and it 

QQ can be used to generate a sampling distribution of sensitivities for a given source 

population. In effect, we derive error bars for the sky-averaged sensitivity curve, 
£f*^ which provide a stringent statistical interpretation for previously unqualified 

• statements about sky-averaged SNRs. As a worked-out example, we consider 

y~ "H isotropic and Galactic-disk populations of monochromatic sources, as observed 

C ) with the "classic LISA" configuration. We confirm that the (standard) inverse- 

£N| rms average sensitivity for the isotropic population remains the same whether or 

r- H not the LISA orbits are included in the computation. However, detector motion 

tightens the distribution of sensitivities, so for 50% of sources the sensitivity is 
within 30% of its average. For the Galactic-disk population, the average and the 
distribution of the sensitivity for a moving detector turn out to be similar to the 
rS isotropic case. 

S3 



PACS numbers: 04.80.Nn, 95.55.Ym 
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1. Introduction 

The new window on the Universe promised by gravitational-wave (GW) astronomy is 
opening, but we cannot as yet see much outside. The eagerly awaited direct detection 
of GWs is likely to be achieved in the next few years with second-generation ground- 
based interferometric detectors such as Advanced LIGO and Virgo [TJ [2J, but at 
the moment we still operate in a predetection era when the planning, design, and 
configuration of future GW observatories is a major research topic. For instance, 
third-generation ground-based instruments such as the Einstein Telescope [3] promise 
much greater reach and yield, but their detailed design is still unsettled. In the space- 
based detection effort, after the ESA-NASA collaboration on LISA [H [5] was halted 
in 2011, ESA has been studying a slightly descoped European-led mission (known 
as NGO, or eLISA (6j [7]) for possible launch at the end of the decade; in late 2011 
the NASA Physics of the Cosmos program issued a Request For Information (RFI, 
[3]) on possible low-cost designs of space-based GW observatories in case NGO is not 
selected, but with a farther implementation horizon. 

The GW sensitivity curve of proposed detector designs describes the GW 
amplitude needed by monochromatic sources to be detectable with a prescribed 
statistical confidence within a given observation time; it is arguably the primary metric 
of the detectors' scientific value, since it figures in the computation of signal-to-noise 
ratios (SNRs, see Sec. [2]), detection rates, and parameter-estimation accuracies, even 
for chirping sources. However, because space-based detectors will move and rotate with 
respect to the long-lived sources that inhabit the low-frequency GW spectrum, there 
is no such thing as a single sensitivity curve — the sensitivity is different for sources 
with different geometries, and even different timing. Due to the motion of the Earth, 
this is true also for ground-based detectors and continuous sources such as deformed 
pulsars. For the sake of convenience, GW scientists have typically resorted to using 
sky-, polarization-, and inclination-averaged sensitivity curves, with the understanding 
that the SNRs quoted for any individual source would be inaccurate, and would usually 
be provided without an error estimate. (While in the GW literature the sensitivity is 
sometimes conflated with the noise spectral density of the detector, these quantities 
have different units and represent different concepts. However, we shall see below that 
the sensitivity can be used to define the effective strain noise of the detector, which 
incorporates the detector's GW response.) 

In this paper we advocate the general approach of computing the geometry- 
dependent sensitivity over a large, physically motivated sample of geometries, and 
then using the resulting distribution of sensitivities to characterize the error, as it 
were, of sky-averaged sensitivity curves. As an example, we work out the "classic 
LISA" sensitivity to monochromatic binary sources that are distributed isotropically 
over the sky or across a realistic Galactic disk. This analysis shows, for instance, that 
for less than 5% of source geometries the sensitivity is worse than twic^f] (or better than 
2/3) the average sensitivity, if we use the "inverse rms" definition of average discussed 
in Sec. [2j and that for 50% of source geometries the sensitivity is within ±30% of the 
average: see Figs. [2] and [4j We also confirm that for an isotropic source population the 
average sensitivity computed, as is standard practice, by neglecting the orbital motion 
of the detector coincides closely with the average sensitivity computed by including 
the motion; for a Galactic-disk population the stationary-detector average sensitivity 

\ Since the sensitivity is defined as the source strength needed to achieve a given SNR, here larger 
is worse, smaller is better. 
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is different, but the moving-detector sensitivity follows closely the isotropic result. 

To enable this analysis, we provide a straightforward, end-to-end recipe 
to compute the exact source-geometry-dependent sensitivity of space-based 
interferometric detectors similar to LISA, including NGO and most of the mission 
concepts submitted to the NASA RFI Indeed, this paper can be considered as 
a useful companion to those studies. The recipe is by no means original: it finds its 
origins in the work of Armstrong, Estabrook, and Tinto [10] , it is analogous to the 
work of other authors |lll I12j . and it expands on the formalism introduced in |13| . 
However, the algorithm described in this paper is self confined and ready to use, and 
it is presented in the context of a careful formal definition of sensitivity. 

This paper is organized as follows. Section[2]is an in-depth discussion of the formal 
and practical definition of sensitivity and of its possible averages; Sections [3] and [3] 
introduce the necessary ingredients for our sensitivity recipe (namely, the noise and 
GW response of LISA-like detectors, including heterodox examples such as "bow-tie" 
and "syzygy" configurations) ; Section [5] puts these ingredients together to compute 
the distribution and average of the classic-LISA sensitivity to monochromatic sources 
(distributed isotropically and across a Galactic disk), and presents sensitivity plots, 
tables, and fits. Last, Appendix A| describes the simple model of the LISA orbits 
used in Sec. [5] while |Appendix B" provides an accurate fit to the average classic-LISA 
sensitivity. 



2. On a general definition of sensitivity 

In signal-processing applications, the detection signal-to-noise ratio (SNR) is defined 
loosely as the strength of a signal divided by the root-mean-square detector noise in 
the same set of data. Intuitively, the SNR is a metric of the confidence, quality, and 
information content of a detection. In the context of a detection scheme, we can make 
precise statements to back up this intuition. 



2.1. Matched filtering 

For instance, in the matched-filtering detection of a signal h immersed in Gaussian 
additive noise n, the SNR of the data s = h + n after filtering by a normalized signal 
template t is given by 

(with the conventions of [14], but our notation), where s(f) and t(f) are the Fourier 
transforms of the detector data and of the template i, respectively; "*" denotes the 
complex conjugate; and S n (f) is the one-sided power spectral density of the noise, 
defined by the ensemble average (n* (f)n(f)) = S n (f)5(f - f')/2 for / > 0. A 
template t is normalized when p(t] t) = 1. 

With this definition, the false-alarm probability of p(s; t) exceeding a chosen 
detection threshold p t hr f or noise alone scales as exp(— /3j hr /2), so higher p(s;t) 
translates into greater detection confidence of a signal h oc t. Indeed, the a priori 
detectability of a signal h is characterized by its optimal SNR \J p(h; h) = SNR opt (/i), 
which is the SNR that would be obtained, on average over noise realizations, using a 
perfectly matching template. 
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2.2. Sensitivity and horizon distance 

These considerations lead to a natural definition of detector sensitivity to a given class 
of signals as the source strength that yields optimal SNR equal to a fiducial detection 
threshold SNR t h r - The definition of source strength is somewhat conventional: in 
matched- filtering searches for GW signals (the concern of this paper), we can write 
h(t) = Aho(t), isolating an amplitude parameter, and then define 

x / A needed for SNR opt 1 SNR thr 
sensu^tyM = | tQ ^ ^ j = (2) 

Because S n appears at the denominator in Eq. ([I]), this definition of sensitivity is in 
effect rms noise divided by rms signal (ho) of fiducial strength. Thus, stronger noise 
makes for worse (numerically higher) sensitivity, which requires stronger signals to 
achieve the same detection quality. Sensitivity can be interpreted readily in terms 
of the horizon distance out to which detection is possible: we write h = (do/d)Aoho 
(with d the luminosity distance and Aq the amplitude at distance do), and define 

horizon(/i ) = . (3) 

For GW applications, a useful quantity to compute is the sensitivity to 
monochromatic sources [with signals h(t) — A cos(27r/o£ + 4>o)], which is related to 
S„(fo) through Eqs. Q and Q by 

sensitivity (/ ) = SNR thr ^^ /o) , (4) 



with T the time of observation. (The equation is dimensionally correct since [S n ] = 
Hz -1 .) This relation highlights two important (if well known) facts: first, the SNR 
for continuous, non-chirping sources grows linearly with the square root of T; second, 
the sensitivity to monochromatic sources contains the same information as S n (f), and 
therefore can be used to reconstruct the SNR of generic sources through Eq. 



2.3. Detector response 

This straightforward idealized picture still requires some elaboration to reflect actual 
GW detectors. This is because the signals recorded by GW detectors are not simply 
proportional to the GW strains at the detector locations; rather, they are junctionals 
of the strains that reflect a variety of physical effects: 

(i) Detector geometry yields the characteristic antenna patterns that describe 
detector response to GWs of given polarization and direction of propagation; 

(ii) Detectors may need calibration — the empirical measurement of the proportional- 
ity and lag of detector output vs. GW strain, as a function of frequency; 

(iii) Detectors may move (with the Earth, or in space) significantly over the timescale 
of GW signals; 

(iv) For wavelengths comparable to (or smaller than) the extent of the detector, 
characteristic time signatures appear in the response. 

This list is not all inclusive. With sufficient generality for our purposes, we are going 
to assume that the detector output Sdet is described (in Fourier space) by 

s d et(/) = R(/):h(/)+n det (/), (5) 
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where the sans-serif quantities denote symmetric trace-free tensors and ":" their inner 
products, R(/) is the detector's response function, and ndot is the instrument noise. In 
the case of ground-based detectors, R(/) includes the detector's calibration function 
[IS] , while the product R(/) : h(/) produces the well-known detector antenna patterns 
F + and F x [TBI Sec. 7.2]. In the case of space-based interferometric detectors, R(/) is 
a complex function that encodes the time-delay structure of the so-called Time-Delay 
Inter ferometry (TDI) observables (see Sec. [3] below), and which is also a function of 
the position and orientation of the spacecraft at the time of observation. 



2.4- Averaged sensitivity 



It is useful to distinguish between the dependence of SNR (and therefore sensitivity) 
on the intrinsic source parameters 6^ (such as binary masses and spins) that determine 
the time evolution and structure of the GW signal, and its dependence on the extrinsic 
parameters f2 M (such as the source's sky position and luminosity distance) that 
describe, as it were, the contingent presentation of the signal to the detector. Now, 
the dependence of sensitivity on the extrinsic parameters is arguably less interestingjf] 
furthermore, we wish to define sensitivity as a function of as few parameters as 
possible, so that we may plot it or fit it easily. Thus, it is often desirable to average 
the sensitivity over the extrinsic parameters. We can, for instance, rms-average the 
optimal SNR over fi, 



<SNR opt (h )) 



0,rms 



\ 



R(/) : ho(/) 



5det(/) 



d/ 



(6) 



and then define the averaged sensitivity as SNR t hr/(SNRo P t(ho))n,rms- Indeed, if the 
GWs can be written as h(/;0, 9) = e(fi) h(f;6), with e(f2) a constant polarization 
tensor (possibly complex), and h(f;0) a function, we can exchange the f2-averaging 
with the integration over frequency, 



(SNRo Pt (h )>o ri 




(7) 



where we have folded the detector's response and Q dependence into the noise, yielding 
the f2-dependent effective strain noise S e g(f,Q) = <Sdct(/)/|R(/) : e | 2 ; and the 
(inverse-rms) averaged effective strain noise S^ s (f), which can be used to define 
a sensitivity to monochromatic sources according to Eq. Q. 

An important example where the factorization h(/;fi,0) = e(O) h(f;6) 
takes place is the dominant quadrupolar emission from inspiraling binaries. For 
monochromatic binaries, for instance, we have h = e + (Q)/i_|_(t) + e x (0)/i x (t) with 

h + = A(l + cos 2 i) cos(27r/ i + (f>o), h x = — 2Acos tsin(27r/ t + (f>o), (8) 
and therefore (for / > 0) 



h(/)= [(l + cos 2 0e++i(-2cos i )e x ] x ^ e **°5(J - f ). 



(9) 



§ Once we detect a source, we will certainly wish to learn about its sky position (e.g., to look for 
electromagnetic counterparts); but we would not normally decide a priori to look preferentially for 
sources from a given sky region. 
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Here / is the GW frequency and A = ^ (^) 5/3 (^Y^ , with d L the luminosity 
distance, and M c = /i 3 ^ 5 Al 2 ^ 5 = (m 1 m 2 ) 3 / 5 /(mi + m.2) 1 / 5 the chirp mass (defined 
in turn by the binary component masses mi and m,2, or by the total mass M and 
reduced mass /i). We should emphasize that this definition of A is conventional; 
another possibility is to take the signal strength as the inclination- (and possibly 
time-) averaged h 2 — h 2 ^ + h\ at the Solar-system barycenter (SSB). This is useful 
because h 2 follows from the binary's energy loss. 

Possible alternative definitions of averaged sensitivity, which do not allow the 
definition of the equivalent of 5™ s (/), include the inverse of mean SNR [i.e., 
SNR thr /(SNR op t)n] and the mean of inverse SNR [i.e., SNR thr x (SNR^)n]. If 
we are using the horizon distance to compute an expected detection rate that is 
proportional to horizon volume, then the relevant average is (SNR 3 , t )J^ 3 (which 

Thorne [17] suggests can be approximated as (3/2)(SNRo pt )J/ 2 for interferometric 
detectors). 

3. Measurement noise of space-based GW interferometers 

Space-based GW observatories similar to LISA @] comprise at least three spacecraft 
orbiting the Sun in formation, with ~ 10 6 km inter-spacecraft distances (for LISA, 
the formation is a quasi-equilateral triangle, with sides ~ 5 x 10 6 km). GW strain 
is measured by monitoring the oscillating distance between freely-falling test masses 
contained by the spacecraft, using heterodyne laser interferometry. In what is known 
as "drag-free" operation, the spacecraft hover around the test masses to protect them 
from external perturbations, performing slight orbital corrections with very precise 
micro-Newton thrusters. The lasers exchanged by the spacecraft are neither split nor 
reflected, as they would be in a laboratory interferometer. Rather, dual one-way laser 
links are established between pairs of spacecraft (i,j), and each link provides a test- 
mass-referenced interferometric measurement of the phase (or equivalently, the 
frequency) of the laser incoming from spacecraft i, compared to spacecraft j's local, 
outgoing laser. Thus, the notations and yji refer to independent measurements 
that compare the same two lasers on different spacecraft (see top left panel of Fig. [I]). 

While the GW response of these measurements is of course central to the 
computation of sensitivity, we shall first discuss their noise, which dictates the 
structure of the TDI observables of LIS A- like observatories ([20]; see [2T] [TS] for 
extensive reviews). For simplicity, in what follows we gloss over several resolvable 
complications that are discussed at length in the literature{{[] we refer our formalism 
to the "classic" non- strapdown LISA architecture^] and we use the simplest possible 
notation (cf. Table I of [22]). 

|| Examples of complications are the cancellation of ultrastable-oscillator (USO) noise, the motions of 
optical benches, the offsets and drifts of laser frequencies, and the Doppler beat notes due to relative 
spacecraft motion. 

1 In the classic architecture, the incoming laser is reflected off the test mass and then compared to 
the local laser on the optical bench; in a strapdown configuration, the incoming laser is compared 
directly to local laser, and a separate interferometric measurement is made of the position of the test 
mass with respect to the optical bench. 
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Figure 1. Left, top: inter-spacecraft measurements yij in a LISA-like three- 
spacecraft configuration. Each arrow represents a laser dj transmitted between 
spacecraft; laser frequencies are compared in the yij measurements, taken at the 
dotted lines. (We do not show the referencing of the measurement to the freely 
falling test masses.) Left, bottom: the unequal-arm Michelson observables of 
time-delay interferometry, centered around spacecraft 1 (observable X), 2 (Y), 
and 3 (Z). TDI observables can be seen as synthesized interferometers, because 
the sequence of time-delayed measurements reproduces the path of light in a table- 
top interferometer, with consecutive measurements (e.g., y2i(t) + j/12 — L21)) 
playing the role of mirrors, and arithmetic subtraction substituting for beam 
splitters and interferometric recombination at photodetectors |18| . Right, top: 
inter-spacecraft measurements yij in a four-spacecraft, eight-link "bow tie" 
configuration. (For clarity, we do not show the local comparison lasers for each 
measurement.) Right, bottom: inter-spacecraft measurements yij in a folded 
"syzygy" configuration |19) . considered under the name SGO lowest in the 2011 
NASA concept study [8] (note however that the and j/32 links would not be 
implemented in that mission concept). 



3.1. Noise in one-way and corner- spacecraft measurements 

The total noise in the individual yij measurements, written in terms of fractional 
frequency fluctuations, is then given by 

= Cij (t - U,(t))-C Tl {t) + y°v(t) - 2yfp(t) : (10) 

for this measurement, which takes place at time t on spacecraft j: 

• represents the frequency noise for the laser sent from spacecraft i to j, which 
is also compared locally with the laser Cji sent from spacecraft j to i; 

• Lij is the one-way transit time of light along the link, which delays the 
instantaneous noise of the incoming laser, as measured in y^; 

• yf™ is the line-of-sight velocity fluctuation of the test mass on spacecraft j; 

• y°J is a sum of optical-path and photon-shot noise for the heterodyne 
measurement. 

"Corner" spacecraft, which participate in laser links with two other spacecraft, include 
two measurement systems, each consisting of a laser, test mass, and optical bench. 
On such a spacecraft (say, j), two additional measurements and Zkj are made to 
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compare the frequency of the two lasers and to remove the noise of the optical fiber 
that connects the two optical benches. These measurements have total noise given by 

zf°™(t) = C k3 (t) - C l3 (t) + 2ylJ{t), (11) 

noi 

Z kj 



z n k f sc (t) = C l3 (t)-C k3 (t) + 2y?™(t). 



3.2. Time-delay inter ferometry 

The yij measurements cannot be used as such to search for GWs, because the laser 
frequency noise of space- worthy lasers would greatly overwhelm all other noise sources, 
as well as the expected response to GWs. Instead, the science data will consist of TDI 
observables: linear combinations of several y and z, appropriately time delayed so 
that all instances of laser frequency noise C cancel out, or rather are suppressed by 
several orders of magnitude. For instance, the unequal-arm Michelson TDI observable 
X (which, as the name suggests, exactly cancels laser noise in configurations with 
constant but unequal arms) is given by 

X(t) = [y 31 (t) + V 31 y 13 (t) + V 13 V 31 y 21 (t) + V 21 V 13 V 31 y 12 {t)] (12) 

- [V2i(t) + V 2l y 12 {t) + V 12 V 21 y 3l {t) + V 31 V 12 V 21 y l3 {t)} 

+ [-V 13 V 31 V 12 V 21 z 31 (t) +V 12 V 21 z 31 (t) +V 13 V 31 z 31 (t) - z 31 (t)]/2 

+ [T>i 2 V 21 V 13 V 31 z 2 i(t) -£> 12 D 21 221 (t) -V 13 V 31 z 2 i (t)+z 2 i(t))] /2, 

where we have introduced time-delay operators T>^ such that for any measurement 
x(t), 

V ljX {t) = x(t- L l3 {t)), VuVijxty^xft-Lijitf-Lkiit-Lijit))^, ... (13) 

Note that in general the delays Lij(t) and Lji(t) in the two directions along the same 
arm are not equal, because of relative spacecraft motion. 

TDI observables can be seen as synthesized interferometers, because their y^ 
components trace the paths and reproduce the phase delays of light across recognizable 
intcrferometric configurations, such as a Michelson interferometer for X — see the 
bottom left panel of Fig. [I] Now, the time-delay operators of Eq. ( 13 ) commute 
only if the armlengths Lij are constant over times comparable to the longest delays. 
If this is not the case, "first-generation" TDI observables such as X cannot cancel 
laser frequency noise to the required accuracy. More complex "second-generation" 
TDI observables have been derived that cancel all laser frequency noise residuals 
proportional to L; however, for the purpose of computing SNRs and sensitivities, 
second-generation observables are equivalent to their first-generation analogs [TSl H2] , 
so we shall consider only the latter in this paper. 

The simplest possible space-based interferometric mission involves two "end- 
mirror" spacecraft exchanging laser links with a single "corner" spacecraft; in that 
case, effectively only one observable equivalent to X can be formed. With the classic 
LISA design, which includes six two-way links between three spacecraft, many different 
TDI observables are possible, some of which use four one-way measurements, some 
five, and some all six. However, all these observables can be reconstructed from a set 
of thre^] linearly independent observables, such as the unequal-arm Michelsons X, 

+ For classic LISA, the accounting of independent observables goes as follows: start with six one- 
way inter-spacecraft measurements yij and three combinations of intra-spacecraft measurements Zij 
(which must be used in pairs to cancel fiber noise), then "solve" for the six lasers that must be 
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Y, and Z centered on the three spacecraft [13]. Two or three such observables can 
be used together to increase SNR and improve GW-parameter determination (this 
insight goes back at least to Cutler [23"]). 



3.3. TDI observables and their noise 



In Ref. [13], one of us (MV) described how the search for TDI observables (and, 
specifically, for a basis {O a } of independent observables) can be turned into a problem 
of linear algebra. This is possible because when the delay operators commute they 
have a convenient representation in the Fourier domain, 



Tpij ■ ■ ■ V u z(t)] = A !3 • • • A kl F[z(t)], with Ai 



(14) 



where T denotes the Fourier transform. For instance, in a four-link-only configuration 
with corner spacecraft labeled by 1, the Fourier-domain dependence of the yij and z^ 
measurements on the laser noises is described by 



2/12 \ 

2/21 
2/13 
2/31 

-*3i)/2 J 



( Ai 2 


-1 








\ 


-1 


A 2 i 
















A13 


-1 










-1 


A 3 i 




V o 


-1 





1 


/ 




(15) 



TDI observables correspond to vectors a of time delays such that a T y = 0, i.e., to the 
solutions of D T a = 0. For the D of Eq. ( 15 ) the only possible solution, except possibly 



for a normalization and overall delay, reproduces the X of Eq. (12 1. 

The non-laser noise spectrum of TDI observables, a key ingredient to building 
sensitivity curves, is obtained by writing one more vector equation that describes 
the y op and y pm content of each y$j and Zij (no delays are involved); multiplying the 
equation by the TDI coefficients a; building the squared complex norm of the resulting 
scalar; and taking its noise-realization ensemble average. Equivalently, we can begin 
with an equation similar to (12 1; replace the yij and with the right-hand sides of 
Eqs. (10 1 and (11); set all the Gy to zero; move to the Fourier domain by switching 
the T>ij to Ay and all noise quantities to their Fourier transforms; and again build the 
squared complex norm and take the ensemble average. Usually (but not necessarily) 
one assumes that all the noises of the same kind can be described as uncorrelated 
Gaussian processes with the same spectral density, 



{frv7 l {f)) = \s op {fW-f)k ]M 



(16) 



(and similarly for y pm ) and that all Ljj = L. For instance, the resulting noise spectrum 
for X is 



S x (f) = 16S op (f)sin 2 x + l6S pni (f)(3 + cos2x)siii 2 



(17) 



where x — 2nfL. For the folded "syzygy" jTH] linear configuration studied as the 2011 
NASA RFI concept SGO lowest [5] see bottom right panel of Fig. [I], the structure 
of X does not change, but the differing armlengths (which we describe by setting 

canceled out. Thus there remain three independent GW-sensitive measurements, although in the 
limit of gravitational wavelengths that are large compared to the armlengths, the symmetrized sum 
of the measurements is insensitive to GWs. 
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L23 = L31 = L32 = L13 = L and L12 = £21 = 2L), result in different trigonometric 
multipliers for the noises: 

&x Zygy (/) = 16 S ° P (/) ( 3 + cos 2x ) sin2 x + 16 SPm (f) ( 5 + 4 cos 2a; + cos Ax ) sin2 ( 18 ) 
where we have also set 5 op = S%f four times larger than all other S??, to account for 
the greater shot noise due to fewer photons available on the longer link. 



3.4. Bases of uncorrelated TDI observables 

In general, the observables in a basis {O a } have correlated noises, since the same 
y°J and y p ™ can appear in different O a s. (Since the Cy have been canceled out, 
we do not need to worry about them.) However, by diagonalizing the matrix S a p — 
(O* (/)O l g(/)) of noise covariances we can find a basis of uncorrelated observables. 
These are especially useful because the total SNR 2 is then the sum of the individual 
SNR 2 . 

For the classic-LISA six-link configuration, three solutions to the equivalent of 
Eq. (151 are given by X and by the Y and Z obtained by cycling indices. The noise 
covariances turn out to be 

S X y = S YZ = S zx = -(8 S op + 32 5 pm ) cos x sin 2 x; (19) 

thus a basis of uncorrelated observables that diagonalize the resulting S a p is given by 
the celebrated A, E, and T, 

A=(Z-X)/V2, E = (X -2Y + Z)/V6, T = {X + Y + Z)/Vs, (20) 

with eigenvalues (i.e., noise spectral densities) 

S a ,e = 8 S" op (2 + cos x) sin 2 x + 16 5 pm (3 + 2 cos x + cos 2x) sin 2 x, (21) 

S T = 16S' op (l-cosa;)sin 2 a; + 128 1 5* pm sin 2 a;sin 4 (a;/2). (22) 

For the less well-known case of a "bow tie" configuration of four spacecraft 
at the vertices of a square, linked by four two-way links, two along diagonals of 
length L, and two along opposite sides of length L/y/2 (see top right panel of 
Fig. [lj, a correlated basis is given (intuitively) by the four unequal-arm Michelsons 
W, X, Y, and Z that can be formed around the vertices 1-4 (defined with signs 
such that W^rX^tY^rZ^tW under the cyclical spacecraft relabeling 
1^2— ^3— s-4— ^1). An uncorrelated basis is given by 

\(W-X + Y-Z), \(W + X~Y + Z), \{W-X-Y + Z), \{W + X + Y + Z), (23) 

with four non-degenerate power spectra, which are somewhat unilluminating functions 
of 5 op , S pm , x, and x/\/2 that we do not show here. 



4. GW response of space-based interferometers 

We turn now to the GW response of LISA-like detectors. The ingredients needed for 
this are the geometry of the constellation (i.e., the individual spacecraft orbits) and 
the time-dependent GW polarization tensors at the locations of the spacecraft. Under 
the very applicable idealization of GWs as plane waves, and under the provably true 
assumption that the constant and orbital-Doppler-induced frequency offsets between 
lasers can be removed once in orbit, the two ingredients come together into the 
Estabrook-Wahlquist two-pulse response for one-way measurements |24) . 

gw U \ _ 1 ' [ h ( p " 1 Z bll Z h fo' *)] ' "'J ioA\ 
Vij [ '~2 l-n-.k ' [ ' 
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where we work in SSB coordinates and where: 

• hij is the unit vector from spacecraft i to j; 

• k is the plane-GW propagation vector; 

• h(f>fc,f) is the TT-gauge GW strain (see, e.g., [16]) at the spacecraft position pt 
and time f; 

• Lij, as before, is the photon time of flight along the arm. 

Thus, an individual GW pulse appears twice in a time series of each y^: once when 
it impinges on spacecraft i, and a time after it has impinged on spacecraft j. The 
GW response of the Zij(t) is effectively zero, because the intra-spacecraft armlength 
is so small. The GW response of the TDI observables is then just a sum of delayed 
y GW (t) terms. 



4-1. Spacecraft orbits 

Since the spacecraft move, the vectors pi and hij and the armlengths L,-j are functions 



of time: specifically, in Eq. (24) they should appear in full as Pi(t— -Ly(t)), Pj(t), hij(t), 
and Lij(t), with the latter two defined implicitly by the light-propagation equation 
Lij(t)hij{t) = Pj(t) —Pi(t — Lij(t)). For many applications, however, it is sufficient to 
evaluate n y -, p iy and at the nominal time of observation t. 

Equation ( p4"| can be used for any type of spacecraft orbits, whether specified 
analytically or numerically. Several of the LISA-likc missions concepts under 
investigation as of late 2011 all adopt variations of the 1 AU heliocentric orbits planned 
for the three LISA spacecraft. These are slightly eccentric Earth- like orbits that are 
inclined with respect to the ecliptic plane (by an ang le ~ y/3/2 x L/AU rad), and that 
lie on three different planes separated by 120°, resulting in a quasi-equilateral-triangle 
spacecraft configuration that is inclined by 60° with respect to the ecliptic. While the 
center of the triangle completes a full revolution around the Sun in the course of a 
year, the triangle itself (as well as the normal to its plane) rotates through 360° in 



a cartwheeling motion (see Fig. Al). The armlengths change slowly, and are always 
equal to one another to ~ 1%. In Appendix A we give simple analytic expressions 



for these orbits; these expressions enforce constant and equal armlengths, and are 
sufficiently accurate for the computation of sensitivity curves. 



4-2. GW source conventions 

Moving on to GW sources, any plane GW incoming along the propagation vector k 
can be described as 

h(i, x) = hssB(i — k ■ x) = e + (k, t{i)h+(t — k ■ x) + e x (k, ip)h x (t — k ■ x), (25) 

where hggB is the value of the metric perturbation at the SSB^jand where h + and 
h x are the two GW polarizations (defined modulo a polarization rotation encoded by 

* The time variable t that appears in h~^~(t) and h x (t) is the time at the origin of the SSB frame, 
which is related to the time in the GW-source frame by a special- and general-rclativistic time dilation; 
it is however expedient to describe the GW using SSB time, including time-dilation factors such as 
cosmological redshift in the parametrization of the GW source. 
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ip), and where e + and e x are polarization tensors such that e + :e + =e x :e x =2, 
that e + : e x = 0, and that e+ • k = e x • k = 0. A possible choice [35] is to take 

/ 1 \ / 1 \ 

e+ = Or ■ -1 • Of, e x = Oj • 1 • Of, (26) 

\ / \ / 

with 

(sin A cos ip — cos A sin (3 sin ip — sin A sin ip — cos A sin j3 cos ip — cos A cos /? \ 

— cos A cos ip — sin A sin j3 sin ip cos A sin ip — sin A sin j3 cos ^ — sin A cos p3 . (27) 

cos j3 sin ip cos /? cos ^ — sin (3 J 

Here (3 and A are the ecliptic latitude and longitude of the GW source, so 

k = — (cos A cos (3, sin A cos (3, sin (3). (28) 

The sky-position parameters /?, A, and ip can be related to the 8, </>, and ip used in Ref. 
[15] and in the Mock LISA Data Challenges HZlUHl 1111301 131] by setting (3 = n/2-6, 
A = <p, and ip = —ip. 

4-3. The GW response to monochromatic signals 

Looking ahead to computing the sensitivity to monochromatic GWs, we can specialize 
Eq. ([24]) to the GW strain of Eqs. ^ and ^ 



GWfn _ i 1 + cos i)h ij -e + -h lj + i(-2cost)n ij -e x -h^ 

Vij 2(1-^4) (29) 

x (^2-Kxf Q {L lj+ k- Pi ) _ e 2 ^/ofc-P J ^e i0o ^(/ - /o) for/ > 0, 

where we have used the Fourier-domain representation of time delays. To get the GW 
response function for a TDI observable, we sum up the individual y^ w (/) for each 
link, multiplying each by the appropriate time-delay coefficient, and take the square 



complex amplitude of the resulting expression. For instance, for the X of Eq. ( 12 1 we 
would write 

A GW (/) = (A 21 + A 21 A 13 A 31 )^ 2 W (/) + (1 + A 12 A 21 ) 2/2 G 1 W (/) (30) 
+ (A 31 + A 31 A 12 A 21 ) 2 / 1 G 3 W (/) + (1 + A 13 A 31 )y 3 G 1 w (/). 



In practice, the difference of two exponentials on the second row of Eq. (29) can 
be factored as e 2 ^°*' w ( e 27ri /°£« - e* ri/ll, ^ ,fi «), and e 2 ^ ifo ^ Pi further decomposed 
by writing pi = R + (pi — R), with R the vector to the "center" of the spacecraft 
formation. The factor e 2m fo k -R then becomes common to all the y GW (/) that appear 
in a TDI observable, and can be dropped (together with e 1 ^ ) when taking the square 
complex amplitude. As a further simplification, in the limit of equal armlengths all 
2-rrfLij = 2irfL = x, and all A t j = e lx . 

Since the spacecraft are always moving along their orbits, frequency-domain 



expressions such as Eq. (29) are only accurate for data segments that are short 
compared to the typical orbital timescale (for LISA-like orbits, a year). However, 
as long as the length of segments is sufficient to resolve the lowest GW frequency of 
interest, the response function can be obtained easily by averaging the individual 
"instantaneous" responses. This approximation neglects the Doppler modulation 
of individual monochromatic sources, which changes results only marginally by 
smoothing out signals across nearby frequencies. 
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For the spacecraft orbits of |Appcndix A| which are described by a simple 
geometric transformation, we can take the alternative, expedient approach of deriving 



the apparent motion of a GW source in the observatory frame [see the end of Appendix 
[A] , and then integrate the stationary-detector sensitivity along that path [see the 
discussion around Eq. (36)]. 



5. Sensitivity to isotropic and Galactic distributions of 
quasi-monochromatic sources 

We now draw together the ingredients introduced in Sees. [3] and [4] to calculate the 
averaged sensitivity of a LISA-like detector to quasi-monochromatic sources, and to 
characterize its variation over sky positions and orientations. Because the detector 
response to GWs is directional, these quantities are only defined with respect to a 
distribution of sources in the sky (i.e., a population): here we consider the simplest, 
isotropic distribution, as well as a realistic Galactic-disk distribution. The latter 
analysis is relevant, for instance, to the sensitivity of LISA-like observatories to 
verification binaries [32] ■ Since verification binaries have known positions, it is possible 
to compute their exact SNR opt (modulo unknown inclinations), but the averaged 
sensitivity and its "error bar" is still useful for quick calculations, or for yet-to-be- 
discovered verification systems. 

Traditionally, the isotropic sky-averaged LISA sensitivity has been computed for 
a stationary detector configuration, reasoning that orbital motion can be neglected 
because, in the LISA "frame," an isotropic source distribution remains invariant as 
the detector's orientation changes along the year, so the sky-averaged sensitivities 
of stationary and moving detectors would coincide. In this section, we compute 
sensitivities in both cases, and show that the argument is indeed correct for the 
averaged sensitivity, but that the resulting sensitivity distributions are rather different. 
For instance, the 5%-quantile sensitivity is ~ 5 times worse than the average sensitivity 
for a stationary detector, but only ~ 2 times worse for a moving detector, whose motion 
effectively smears away the "blind spots" of the stationary detector. 

This section is organized as follows: in Sec. |5.1| we describe the details of our 
calculation, including our assumptions for the instrument noises, for the residual 
confusion noise from the Galactic white-dwarf background, and for our Galactic-disk 
population. In Sec. |5.2| we plot the resulting distributions of sensitivities, discuss their 
features, and provide frequency-dependent fits. 



5.1. Computing the sensitivity 

We now proceed to compute the LISA-like sensitivity to monochromatic-binary GW 
signals for a single TDI observable X. We consider two detector geometries: i) a 
stationary equilateral-triangle constellation with 5 Mkm arms, lying in the ecliptic 
plane and centered at the Sun; ii) a rotating constellation in solar orbit, with the 
LISA-like spacecraft trajectories described in | Appendix A| 

For both cases we randomly select 10 5 source positions in the sky, distributed 
either isotropically (uniformly in sin/3 and A), or across a thin-disk Galactic model of 
density 

dN 

oc r exp(— r/2.5kpc) sech(z/300 pc) (31) 
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( S op)l/2 ) physical (SP" 1 ) 1 / 2 , physical 

S op , fractional frequency S pul , fractional frequency 

Mock LISA data challenges 
20 x 10" 12 mHz- 1 / 2 3 x 10~ 15 [1 + (^)' 2 } 1/2 

1.8 x lO-^/^Hz- 1 2.5 x io- 4 8 [1 + ( Jh 2 _ ) -2 ]/ -| Hz -i 

LISA science requirements (2007-2011) 
18 x 10- 12 [1 + )- 4 ] 1 / 2 mHz-V2 3 x Hr 15 [l + ( I f^)- 1 ] 1 / 2 [l + ) 4 ] 1/2 ^72 

1.42 X 10- 37 [1 + (^f^)- 4 ] / 2 Z Hz- 1 2.53 X l(r«[l + (^"^[l + (^^) 4 ] /£? Hz" 1 

ESA NGO (2011, final) 
(5.25 X 10" 23 + 6.28 X lO" 23 ) 1 / 2 m Hz" 1 / 2 4.6 X 1Q- 15 [1 + (t^)- 1 ] 1 / 2 — " 



m s 



(2.31 X 10" 38 + 2.76 X 10- 38 )/ 2 z Hz" 1 6 X 1Q- 48 [1 + ( jfe)- 1 ]/^ 2 Hz" 1 



iO- 4; J H Z !/2 
H 

10- 



Table 1. Optical-path and proof-mass noise components for three models of 
LISA-like GW observatories. Noises are given in physical units and in the 
fractional-frequency fluctuations units used in this paper. Note that for ESA 
NGO, the reference armlength is 1 Mkm rather than LISA'S 5 Mkm. The two 
S op components given for ESA NGO correspond to shot noise (which scales with 
armlength, laser power, and telescope size) and to other fixed optical noises. 



[33] . where r and z are the Galactic cylindrical radius and height (they can be 
converted to ecliptic j3 and A following, e.g., [M]). For each source position we also 
randomly select the inclination 1 and polarization angle ip from isotropic distributions. 

We use the instrument noise of Eq. (17), with the assumptions on optical-path 
and proof-mass noises given in the 2011 LISA requirements [35], and spelled out in 
Table [T] where they are compared with the parameters used for the Mock LISA Data 
Challenges [27J EM CM IH1 EE] and for the ESA 2011 "NGO" study 0[7J. 

The instrument, however, is not the only source of noise that affects the sensitivity. 
LISA-likc detectors with sufficiently long baselines and low instrument noise will be 
subject to an unresolvable GW confusion foreground from white-dwarf binaries in the 
Galaxy (see |36| for a recent treatment and extensive references). In this paper we 
adopt the confusion- foreground model developed in [36] . 

'/3dA x r ° 1«~2„ 5 „2. 



^° nf (/) = 1-4 x 10~ 45 (//Hz)- 8 / 3 Hz" 



x tanh° 



2 d/ 



x 



20 (1 + 0.6cc 2 



(32) 



(here expressed as fractional frequency fluctuations for the TDI observable X), where 
a~ 1, j3 ~ 1/yr, dN/df = 5 x lO^Z/IiVr 11 / 3 Hz" 1 , and x = 2tt/L. The first term 
in Eq. (32) is a fit of the GW-strain power spectral density of unresolved Galactic- 
binary sources, in the mid-frequency regime where systems are evolving because of 
GW emission alone; the second term is a fit of the transition to the higher frequencies 
(above a few mHz) where binary signals are sufficiently spaced apart in frequency 
domain to be resolved and subtracted individually; the third term is a sky-averaged, 
approximated transfer function from equivalent strain to X. The total X noise is then 
given by Sx(f) + Sx nl (f), and is comparable to (but measurably different from) the 
often-seen expression proposed by Barack and Cutler [37] . 

When we compute moving-detector sensitivities, we include the strong seasonal 
variation in the amplitude of the confusion foreground caused by the evolving 
orientation of the LISA-like formation with respect to the Galaxy. Following Ref. 
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[36] , and consistently with the orbits of Appendix A| we set 

S°? nf (/, t) = S%°* (/) x [1 - 0.21 cos(27rt/yr) - 0.45 cos(47rt/yr) (33) 

- 0.05 sin(27rf/yr) + 0.01 sin(47rf/yr)] . 

For each source geometry (sky position, inclination, and polarization, denoted 
collectively by f2), we compute a full sensitivity curve from / = 10~ 4 to 10" 1 Hz. 
Specifically, we insert the binary quadrupole-emission polarization tensor [ the first 
bracket of Eq. ^] into the frequency-dependent GW response X GW of Eqs. (29) and 
( 30 ) , using the e + and e x of Eqs. (261 and (27). The effective strain noise (i.e., the 
TDI noise over the Sl-dcpendent squared TDI response function) is then 

c (f.Q\- S *(f) + S x ni (f) (OA'S 

b eS {f,ll)- i X GW(y ; o)|2 ' W 
and the corresponding f2-dependent monochromatic-signal sensitivity is 

scnsitivity(/ ; fi) = SNRthr V^pM , (35) 



In the next section, we study the distribution of this sensitivity function over our 
sample of geometries f2. We define the average sensitivity function from the inverse- 
rms averaged effective noise 5™(/) = (S e g(f; Q)^ 1 )^ 1 . We always set SNRthr = 5 
and T = 1 yr. 

As mentioned at the end of Sec. |4.3| and in |Appcndix A[ we handle the moving- 
detector case by imposing an apparent motion tt'(t;Q) on the sources, and integrate 
the resulting instantaneous effective strain noise, 



' S x (f) + ST l (f,t) 

\x^(f-,n'(f,n))\ 2 



s „ {f . n)= dt ~x w , / (36) 



In practice, we compute this integral as a discrete sum of terms evaluated at N pos 
positions along the year-long LISA-like orbits. We find that N pos = 26 provides a 
sufficiently accurate approximation, with a typical error -C 1% (a larger A^os will be 
required at frequencies higher than we consider here). 



5.2. Sensitivity curves: distributions and averages 

Figure [2] shows the distribution and average of the sensitivity, computed over the 
isotropic Q sample as discussed in the previous section, for the stationary (top panel) 
and moving (bottom panel) detector configurations. At any frequency, the vertical 
gray-level profile encodes the quantiles of the sensitivity distribution: for instance, 
the sensitivity would lie above the curve labeled "0.020" (or below the curve labeled 
"0.980") for 2% of source geometries. The continuous and dashed light-colored lines 
mark the (inverse-rms) average and median sensitivities, while the histograms on the 
right show the distributions at three representative frequencies, for which we also give 
numerical values in Table [21 

For the stationary detector, the distributions have long tails: the extremal values 
in our O sample can be as much as three orders of magnitude worse than the 
average, corresponding to "blind spots" in the sky, caused by unfavorable orientations 
of the spacecraft constellation with respect to the polarization of the incoming 
GWs. However, these tails are thin, as shown by the position of the 2% and 98% 
quantile curves. Over the entire frequency range, the average sensitivity corresponds 
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log 10 (f/Hz) p(sensitivity) 



Figure 2. Top panel: Distribution of sensitivities for a stationary LISA- 
like detector, over a sample of 10 5 source geometries selected isotropically across 
the sky. Gray levels map the quantiles of the distribution (i.e., the fraction of 
the total sample for which sensitivity is better than plotted); the four black 
curves correspond to 2%, 10%, 90%, and 98% quantiles. The inverse-rms 
averaged and median sensitivities are plotted by the solid and dashed yellow 
lines, respectively. The panel on the right shows histograms of sensitivities for 
representative frequencies of 0.3, 10, and 100 mHz; these are fit very well by 
the skew-normal distributions of Eq. |37| , plotted by the red curves. Bottom 
panel: Distributions of sensitivities for a moving LISA-like detector. The 
histograms in the right pan el a re fit well by the modified, double-peaked skew- 
normal distributions of Eq. (| 39 [> - 
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sensitivity, stationary detector (absolute/relative to average) 





0.3 mHz 


10 mHz 


100 mHz 


5% 


5.26 x 10" 


-21 


/ 4.88 


2.58 x 10" 


-23 


/ 4.65 


1.60 x 10" 


-22 


/ 4.71 


10% 


3.52 x 10" 


-21 


/ 3.26 


1.78 x 10- 


-23 


/ 3.20 


1.14 x 10" 


-22 


/ 3.37 


average 


1.08 x 10" 


-21 


/I 


5.54 x 10" 


-24 


/I 


3.40 x 10" 


-23 


/I 


median 


1.32 x 10" 


-21 


/ 1.22 


6.82 x 10" 


-24 


/ 1.23 


4.46 x 10" 


-23 


/ 1.31 


90% 


6.96 x 10" 


-22 


/ 0.65 


3.57 x 10- 


-24 


/ 0.64 


2.18 x 10" 


-23 


/ 0.64 


95% 


6.01 x 10" 


-22 


/ 0.56 


3.08 x 10" 


-24 


/ 0.56 


1.84 x 10" 


-23 


/ 0.54 



sensitivity, moving detector (absolute/relative to average) 





0.3 mHz 


10 mHz 


100 mHz 


5% 


2.00 x 10" 


-21 


/ 1.93 


1.06 x 10" 


-23 


/ 1.91 


6.84 x 10" 


-23 


/ 2.01 


10% 


1.82 x 10" 


-21 


/ 1.76 


9.56 x 10" 


-24 


/ 1.72 


6.29 x 10" 


-23 


/ 1.85 


average 


1.04 x 10" 


-21 


/I 


5.55 x 10" 


-24 


/I 


3.40 x 10" 


-23 


/I 


median 


1.16 x 10" 


-21 


/ 1.12 


6.21 x 10" 


-24 


/ 1.12 


3.87 x 10" 


-23 


/ 1.14 


90% 


7.29 x 10" 


-22 


/ 0.70 


3.90 x 10" 


-24 


/ 0.70 


2.39 x 10" 


-23 


/ 0.70 


95% 


6.85 x 10" 


-22 


/ 0.66 


3.68 x 10" 


-24 


/ 0.66 


2.16 x 10" 


-23 


/ 0.64 



Table 2. Sensitivities of an isotropic population of monochromatic 
sources for stationary and moving LISA configurations listed as (absolute 
scnsitivity)/(sensitivity relative to average) for the 5%, 10%, 90%, 95% quantiles, 
the average, and the median sensitivities for source frequencies of 0.3, 10, and 100 
mHz. 



approximately (±5%) to the 67% quantile curve — thus, about 33% of all sampled 
source geometries have better (lower) sensitivity. 

For the moving detector, the distributions are tighter, with "error bars" a factor 
~ 2 smaller, and almost no tails: since the orientation of the detector with respect 
to any source changes across the year, any blind spot is smeared out. The average 
sensitivity corresponds approximately to the 60% quantile curve. As expected, the 
stationary- and moving-detector average sensitivities coincide very accurately. In 
Appendix B we provide a fit that is therefore valid for both. 

We find that the stationary-detector sensitivity distribution is fit well, at any 
frequency in the range shown, by the skew-normal distribution 



Skew(z; A, £, a, oS) = 



A 



exp 



1 + erf 



V2 



where z = (y — and where 
JJ = log'i 



sensitivity(/; f2) 



(37) 



(38) 



' 10 (sensitivity(/; n)) invcrscrms ' 

The parameters of Eq. p7] ) specify the shape (a), scale (w), center (£), and 
normalization (A) of the skew normal. These are all functions of frequency and 
are shown by the solid curves in Fig. [3j For [0.1, 10] mHz, the parameters are 
approximately constant, with a ss 2.5, u ~ 0.38, £ ~ —0.19, and A w 0.98. 

The moving-detector sensitivity distribution is qualitatively different than 
the stationary-detector case, and has a somewhat peculiar double-peaked, nearly 



Non-sky-averaged sensitivity curves for space-based GW observatories 18 



6 
4 
2 


0.65 
0.60 
0.55 
0.50 
0.45 
0.40 
0.35 




0.30 
10 




it 
o 



10 3 10' 
/[Hz] 



10" 





1.0 
0.5 
0.0 
-0.5 
-1.0 



-1.5 
10 



















! 









10 J 10' 
/[Hz] 



10" J 



Figure 3. Skew-normal and modified skew-normal fit parameters a, oj, £, A, 
bo, and b2, as functions of source frequency /, for stationary (solid) and moving 
(dashed) LISA-like detectors and an isotropic monochromatic-source population. 



1 



(39) 



symmetrical shape. Surprisingly, it is fit rather well by a slight variation of (37) 
that replaces the argument of the exponential with a fourth-order polynomial, 

™ , h , x A ( (z 2 - b 2 z - b ) 2 
bkew(z;A,£,a,u},b ,b 2 ) = — 7= exp 

with the two additional structure parameters b and b 2 - All six moving-detector fit 
parameters are shown as functions of frequency by the dashed curves of Fig. [3| Unlike 
the stationary-detector case, the parameters are sensitive to the presence of confusion 
noise from the Galactic white-dwarf foreground, as shown by the characteristic bumps 
below a few mHz. 

Last, we consider a sample of source geometries distributed according to the 
Galactic-disk population of Eq. (31). As mentioned above, this analysis is relevant to 



the sensitivity of LISA-like detectors to verification binaries, which are in fact rather 
good realizations of monochromatic GW sources. Again, Fig. [4] shows the distribution 
and average of the sensitivity, for the stationary (top panel) and moving detector 
(bottom panel) configurations. For the stationary case, we see that the deviations from 
the average sensitivity can be much larger than for the isotropic population of Fig. 
[2] The average sensitivities themselves are rather different: compare the solid yellow 
line with the dotted black line. For a moving detector, sensitivities are concentrated 
much more tightly around the average, and again exhibit double peaks. Furthermore, 
the average sensitivity is essentially indistinguishable from the isotropic case. 

We were unable to find high-accuracy fits to these distributions. The skew normal 
of Eq. (37) works well for a stationary detector at frequencies less than 10 or 20 mHz, 
but fails dramatically at higher frequencies. The modified skew normal of Eq. (39) 
works reasonably well for a moving detector at all frequencies, although it consistently 
fails to resolve accurately the sharp peaks of the histograms seen in the bottom-panel 
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Figure 4. Top panel: Distribution of sensitivities for a stationary LISA- 
like detector, over a sample of 10 5 geometries for sources distributed across the 
Galactic disk of Eq. \3l\ . All curves are to be interpreted as in Fig. [2] the 
additional dotted black curve plots the average sensitivity for the isotropic source 
distribution (the yellow solid curve of Fig. |2j. Bottom panel: Distribution of 
sensitivities for a moving LISA-like detector. 
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Figure 5. Modified skew-normal fit parameters a, lo, £, A, f>o> an< i ^2, as 
functions of source frequency /, for moving LISA-like detectors and a population 
of monochromatic sources from the Galactic disk. 



of Fig. |4j Nevertheless, we do recommend that Eq. (39 1 be used to provide error bars 
for verification-binary sensitivities. In practice, this can be attained by looking up the 
values of the modified skew-normal parameters at the desired source frequency in Fig. 
[5] and substituting those values into Eq. ( 39 ) . This yields a function of z — (y — £) /uj 
where y is given in Eq. ( 38 ) and the inverse-rms sensitivity at the desired source 



frequency can be found for the appropriate spacecraft configuration from Figs. [2] and 
[4] or from our fitting functions in Eq. (B.2|. If using the latter one must add the 
confusion noise to the function defined in Eq. (B.2). The resulting function of the 



sensitivity yields the approximate distribution at the given frequency. 



6. Summary and applications 

In this paper we present a straightforward end-to-end recipe to compute the non- 
sky-averaged sensitivity of interferometric space-based GW detectors. The recipe 
incorporates the motion of the spacecraft and the seasonal variations in the partially 
subtracted confusion foreground from Galactic binaries. By evaluating the sensitivity 
over a sample of source geometries, we are able to provide a stringent statistical 
interpretation for previously unqualified statements about sky-averaged SNRs. As an 
example, we consider the "classic LISA" sensitivity to monochromatic GW sources 
that are distributed isotropically or according to a Galactic-disk population. 

For an isotropic source population, we confirm that the inverse-rms sky-averaged 
sensitivity is the same whether or not the motion of the detector is taken into 
account. However, the stationary-detector calculation overestimates the variance of 
the sensitivity distribution (see Fig. [2] and Table [2]), because with orbital motion the 
"blind spots" of the detector trace a path across the sky, and are averaged out over 
a year. For a Galactic-disk population, the stationary- and moving-detector average 
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sensitivities are rather different (see Fig. |4j) ; indeed, the former will depend strongly 
on the arbitrary orientation of the stationary detector. However, the moving-detector 
curve coincides closely with the isotropic average sensitivity. Thus, the standard 
approximate sensitivity expression described in Appendix B (where it is also improved 
with the correct high-frequency behavior) can be used for both populations. The 
modified skew-normal fit developed in Sec. |5.2| can be used to describe the statistical 
distribution of sensitivities as a function of frequency, thus quantifying the errors in 
future studies of GW-source detectability. Approximate fit parameters can be looked 
up m Fig. [3] or Fig. [5] and substituted in Eq. ([39]). 

Our recipe can be easily adapted, and our analysis repeated, for different detector 
geometries and orbits; indeed, we have done so for the mission concepts submitted 
to the NASA RFI. For eLISA/NGO, we have confirmed that the distributions of 
sensitivities are fit well by Eqs. [37] and [39] albeit with different fitting parameters. 
Our model of NGO included a constellation with the orbits of | Appendix A| but with 
1-Mkm armlengths, the instrument noises listed in Table [l] plus confusion noise given 
by Eq. @ with a ~ 0.7, (3 ~ 1.2/yr [H]. 

Sensitivity curves can also be obtained for different types of sources, such as 
binary black-hole inspirals. However, since these signals are broadband, it is more 
appropriate to derive the horizon distance (e.g., as a function of total or chirp mass) 
rather than the sensitivity as a function of frequency. Furthermore, the relative timing 
of the inspiral and detector orbits also enters the calculation, and needs to be included 
in the statistical sample. To facilitate such studies, we are providing a Mathematica 
notebook that implements our calculation, which can be easily modified for different 
detectors and sources, which includes fitting parameters for "our" LISA and NGO, 
and which can be found in the CQG supplementary materials for this paper, and at 



http : //www . vallis . org/publications/ sensitivity 
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Appendix A. A simple geometric model for the motion of LISA-like GW 
interferometers 

We begin by placing the spacecraft along the vertices of an equilateral triangle, working 
in a moving, rotating coordinate frame where the spacecraft are at rest |25j : 

p rcst = (£/ 1 /3)(_ CO s2o- i ,sin2cr i ,0) with ^ = 3tt/2 - 2(i - 1)tt/3, (A.l) 

where L is the common armlength. Next, we switch to an inertial, ecliptic coordinate 
system, with origin at the Solar-system barycenter (SSB). In this frame, the center R 
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^moving 



Figure Al. Left: Orbital motion of a LISA-like GW interferometer, as described 
by the simple geometric model of ]Appendix A| shown in an ecliptic SSB coordinate 
system. Right: Apparent motion in the rotating detector frame of sources at 
ecliptic sky coordinates (/3 truo , A truc ). 



of the LISA-like formation moves in the x— y plane on an SSB-centric circular orbit 
with radius 1 AU, and the three spacecraft rotate around R: 

Mt)=R + 2 - P r\ (A.2) 
where R = (1 AU) x (cos 77, sin 77, 0), and O2 is the rotation matrix 

(sin r] cos £ — cos 77 sin £ sin £ — sin 77 sin £ — cos r\ sin £ cos £ — cos 77 cos £ \ 
— cos i] cos £ — sin rj sin £ sin £ cos 77 sin £ — sin 77 sin £ cos £ —sin 77 cos £ .(A. 3) 
cos C sin £ cos £ cos £ — sin £ / 

The functions 7/ = fit + 770 and £ = — fit + £0 that appear in this matrix represent 
the true anomaly of the spacecraft formation and the relative phase of the spacecraft, 
respectively; £ = —ir/6 sets the 60° inclination of the triangle with respect to the 
ecliptic. The resulting orbits are shown in the left panel of Fig. |A1| 

For simplicity, in this paper we set 770 = £0 = 0; however, to map these orbits 
to those used in the LISA Simulator and in the Mock LISA Data Challenges 
[27J US 123 EM EI] , we need to relabel spacecraft 1, 2, and 3 as 0, 2, and 1, and then 
set 770 = K, £0 = 37r/2 — k + A, where n and A are the parameters defined below Eqs. 
(56) and (57) of Ref. [26]. 

The apparent motion in the spacecraft rest frame of a GW source that is located 
at ecliptic latitude (3 and longitude A, and that has principal polarization axes oriented 



along ip (see Sec. 4.2 1 can be derived by requiring that the geometric products that 



appear in Eq. (24 1 remain unchanged as the spacecraft move along their orbits. 



Specifically we require that 

.moving . £true = (0 2 fi- st ) . fct™ = 7^ cst • (O^feff 6 ) = n\f ■ fc™ 8 ; (A.4) 

by equating the components of ^™ ovlng = _ ( CO s A' cos f3' , sin A' cos j3' , sin p' ) with those 

of 0% k\j ne = —O2 ■ (cos A cos (3, sin A cos /3, sin (3) we can then solve for (3'(t) and X'(t) 
as functions of /?, A, r](t), and £(rj) (more precisely, we use the three equations to solve 
for cos A', sin A', and sin/3'). Some examples of the resulting apparent motion are 



shown in the right panel of Fig. Al To obtain the apparent variation ip'{t) of the 



polarization angle we need an additional equation that we can obtain from 

.moving , [ f-e + (O t 1 ruc ) T ] • 77™ s = h]f ■ [Or^e+iOT^f] ■ nlf, (A.5) 
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Figure Bl. Top: inverse-rms averaged sensitivity for an isotropic distribution of 
source ge omet ries for the "classic LISA" configuration (black, solid), the standard 
fit of Eq. ( |B . 1 [> (slightly improved by rescaling S° p ) (black, dashed), and the better 
fit of Eq. \B.2\ (red, solid). As usual, here we set SNR t h r = 5, T = 1 yr. Bottom: 
Fractional error of the improved fit. 

which implies 

Or vins (/?', A', V') = ■ 0\ est (P,\,ip). (A.6) 

This is a matrix equation, but any two components that depend on cos ip' and simp 1 
can be used to solve for ip'(/3,\,i/j) after replacing the functions found above for f3' 
and A'. 



Appendix B. A fit to the averaged sensitivity for stationary and moving 
LISA-like detectors 



The LISA science requirements document [35] provides the standard approximation 
for the inverse-rms averaged effective-strain noise for LISA-like detectors, 



S 



hid 



cir 



(/) 



sin^ 60° 



1 



2fL_ 
0.41c 



4 5'P m (/) 2 



+ S op (f) 2 \/L 



(B.l) 



where 60° is the angle between the interferometer arms, c is the speed of light, L is the 
interferometer armlength (5 x 10 6 km for LISA), and the optical-path and proof-mass 
noise components are given in Table [I] (for this expression, one should use their values 
in physical units). Equation ( |B.l I does not include the Galactic confusion foreground. 
The sensitivity can then be obtained as SNR thr y/S^/T. 



Equation (B.l I does not model the oscillatory response- function behavior evident 
at frequencies larger than about c/2L (30 mHz for LISA). We can do better by 



replacing S op with 0.964 S op , and by multiplying Eq. (B.l ) by the correction function 



S. 



rms,fit 
eff 



rri'ms,std 

— ^eff X 




Skew{f;A,t,a,u) 
2/ a 2 (/ - a 3 )) + a 4 



/ < 0.015 Hz 

0.015 < / < 0.0375 Hz, (B.2) 
/ > 0.0375 Hz 
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where a = 0.9587, a_ = 0.3319, a 2 = 104.7, a 3 = 0.01501, and a A = 0.8309. 
function Skew(/; A, £, a, co) is the skew-normal distribution 



The 



Skew(f;A,£,a,u) 



.4 



'2it 



exp 



U-if 

2a; 2 



1 + erf 



V2lo 



(B.3) 



with A = 0.009752, a = -1.3908, £ = 0.03093, and ui = 0.006365. Figure [bT] shows 
the sensitivity curves resulting from this fit and from Eq. (B.l), compared with the 
numerically derived inverse rms sensitivity (the yellow solid curve of Fig. [2] and of 
the bottom panel of Fig. EJ). The fit is accurate to better than 5%: the skew-normal 
distribution fits the first "hump" between 15 and 37.5 mHz to 1%, whereas the sinusoid 
fits the remaining oscillations to 3%. 
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